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Abstract 

Effects of disorder on electron-doped iron pnictides are investigated systematically based on self-consistent 
Bogoliubov-de Gennes equations. Multiply impurities 'with same scattering potential (SP) are randomly distributed 
in a square lattice. Probability distribution functions of normalized order parameters for different impurity concen¬ 
trations 5imp, different electron doping concentrations 6 are investigated for given SPs. Samples are found to be very 
robust against weak SP, in which order parameters do not have qualitative change even at very large 6imp- While 
strong SP is able to easily break down the order parameters. For moderate SP, variations of order parameters on 
and around impurities strongly depend on d, however the distribution functions of normalized order parameters have 
similar behavior as 6imp increases. Compared with superconducting (SC) order, the magnetic order is more sensitive 
to multi-impurity effect. The spatial spin density 'wave pattern has already been destroyed before the system loses its 
superconductivity. Dependence of SC order on temperature is similar to that of impurity-free case, with the critical 
temperature being remarkably suppressed for high Simp- 
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1. introduction 

The discovery of the new family of layered iron- 
based SC materials has attracted much attention since 
their discovery. (11] The parent compounds of iron pnic¬ 
tides are bad metal at low temperatures and exhibit a 
spin-density-wave (SDW) order. Upon doping either 
electrons or holes into the system, SDW will be sup¬ 
pressed and superconductivity emerges. Disorder thus 
cannot be avoided in the crystal gro'wth procedure 'with 
intentional doping. By substituting the iron ions with 
transition-metal (TM) elements Co orNi, (^.Sfllll] su¬ 
perconductivity can be realized and the dopant ions act 
as impurities at the same time. 

Anderson’s theorem tells us that the i-wave super¬ 
conductivity is insensitive to small SP impurity. While 
in cuprates, impurity effect on the local density of states 
has been used as a prominent feature of c/-wave super¬ 
conductor. 10,13,0] For iron-based superconductors, im- 
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purity scatte ring effect has been theoretically studied in¬ 
tensively ( 9, 13 H. 12 12 i 3 15], but mainly by con¬ 
sidering a single impurity. 

In TM doped cases, doping density is expected to 
be equal to im puri ty concentration. X-ray absorp¬ 
tion experiment lllbll and angle-resolved photoemission 
spectroscopy experiment jlTI] show that TM substitu¬ 
tion effect as well as multi-impurity effect are rich in 
physics. Experiments have measured the variation of 
superfluid density as function of temperature in crys¬ 
tal of Ba(Fei_xCox) 2 As 2 and in film Ba(Fei_xCOx) 2 As 2 . 
The single crystal lIlSll sho'ws a BCS-like behavior, 
■while the epitaxial film lll9ll has a wide range of linear 
dependance over temperature. The difference bet'ween 
the crystal and film samples is usually considered that 
films are dirtier with large disorder. 

Impurity effects are sensitive to SC order symmetry. 
One commonly believed idea is that the SC order is 
formed via spin fluctuation, leading to i±-wave iaHi 
22 . 3l pairing symmetry in iron pnictide. In this paper 


we address the robustness of i±-wave superconductors 
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against multi-impurity and only concern electron doped 
cases. 

By employing a phenomenological model based on 
Bogoliubov-de Gennes (BdG) equations, we study sys¬ 
tematically the multi-impurity effect in iron pnictides. 
The impurity Hamiltonian is described by an on-site 
point-like (5-function potential, with SP strength Ump 
ranging from weak to strong. Multiple impurities have 
interference effect that can smear out the effect of single 
impurity. Although order parameters may be suppressed 
or enhanced on and around impurities, average values 
of order parameters decrease as 6 imp increases. And as 
a whole lattice they obey a statistical distribution func¬ 
tion. We will study the probability distribution func¬ 
tions of normalized order parameters for different Simp, 
different S and different Uimp- With increasing Simp, the 
height of the peak of distribution functions falls down 
at the beginning and then goes up for more larger Simp 
corresponding to the vanished order. 

For weak SP, average values of the order parameters 
do not have qualitative change from that of the impurity- 
free case even when Simp is very large, suggesting that 
the system is very robust against weak multi-impurity. 
For moderate SP, SC order has finite values at Simp - 
0.2, however, the stripe-like magnetic structure is de¬ 
stroyed. Strong SP breaks down the order parameters, 
for Uimp = 5.0, the SOW pattern immediately disappear 
at Simp - 0.0125. 

The rest of the paper is organized as follows. In 
Sec.m we introduce the model. In Sec. |2 we inves¬ 
tigate the small SP cases, and in Sec. |4] and Sec.|5]we 
discuss the moderate and strong SP cases respectively. 
Finally, we provide a summary and conclusion. 


2. model and formalism 


We use a two-orbital four-band phenomenological 
model lUlllilllllill^lIilliSto describe the 


Hamiltonian of iron pnictide. This model takes the dxz 
and dyz orbitals of the Fe ions into account and each unit 
cell accommodates two inequivalent Fe ions. It repro¬ 
duces qualitatively the energy band structure observed 


eiy me energy bai 

by experiments ^ 33,3 ^l36l]. T he obtained re¬ 
sults are consistent with the ARPES 03711 . neutron scat¬ 
tering lli^l experiments. 

The total Hamiltonian is defined as // = Hbcs +Hi„t + 
Himp. Here Hbcs is the BCS-like Hamiltonian, which 
includes the hopping term and the pairing term. The in¬ 
teraction Hamiltonian includes the on-site Coulomb 
interaction U and Hund’s coupling Jh at the mean-field 


level. They are described by 124, 25,26, 13 28, 29 


Hbcs - ^ (diy'jvcly,^Cjva- + h.C.) p ^ c!^,^Cipcr 

iv'jvfT iv'cr 

+ 2 (Ai,-j,4,^c.'„-+/r.c.) (la) 

iv'jvcr 

Hint — H ^ ^ {f^iva}^ivcr + — 3Jh^ ^ ^ (^iy'o-)^ivo- 

i,v,cr^(T 

+ (U-2Jh) ^ <niv'o->«iy<r- (lb) 

i,v'^v,cr^d' 


In the above equations tiiya- is the electron density oper¬ 
ator at site i, orbital v, and spin cr, and p is the chem¬ 
ical potential. ti represents the nearest-neighbor (nn) 
hopping between the same orbitals on Fe ions, 12 ( 13 ) de¬ 
notes the next-nearest-neighbor (nnn) hopping between 
the same orbitals mediated by the up (down) As ions, 
f 4 is the nnn hopping between different orbitals. See a 
schematic illustration in j^ . We adopted the hopping 
parameters as fi, 2 , 3,4 = 1,0.4,-2.0,0.04, respectively. 
The energy and SP are measured in unit of f i. The pair¬ 
ing order Ay = = ^{ciy'^c-^yi - QvqCjvr) is chosen 

as the nnn intra-orbital pairing with v' = v and Vsc - 
1.5. This kind of pairing is consistent with the s±-wave 
superconductivity [23 HI 39, ^ 411 and has been 


widely used in previous studies EIM12711281 
Impurity Hamiltonian Himp is given by 


liH. 


Himp ~ 






( 2 ) 


where m is the impurity site acting as a scattering cen¬ 
ter with Uimp varying from weak to strong. Our numer¬ 
ical calculation is performed on a 40 x 40 square lat¬ 
tice with the periodic boundary conditions. For fixed 
Simp, S and Uimp we self-consistently solve the BdG 
equation iffsi E3] to obtain SC order Ai = 
and magnetic order mi = 2 2v(”ivT “ ”ivi), where t is 
the four nnn diagonal links from site i. In the absence 
of impurities, m, shows a well-known stripe-like anti- 
fetTomagnetic structure with two nearest rows having 
opposite nii. 

In the calculations, the positions of the multiple im¬ 
purities are spatially random. In order to obtain reli¬ 
able probability distribution functions, for a fixed con¬ 
centration Simp, we carry out several numerical calcu¬ 
lations of different configurations. The probability dis¬ 
tribution functions of SC order pA and magnetic order 
Pl„,l are then statistically analyzed. pA(P|m|) is denoted as 
the probability of Ai/Ao(|OTi|/|mo|), and here we define 
Ao and |mo| as the order parameters in the correspond¬ 
ing impurity-free case. Summation of pa and p|„,| are 
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Figure 1: (color online) Plots of pA and p\,„\ for Uimp = 0.5 and Simp = 
0.025,0.1 . Panels (a),(b),(c) coiTespond to d = 0.04,0.08,0.16, re¬ 
spectively. The X coordinate of the black(blue) solid square in each 
panel denotes the average value of A/ for Simp = 0.025(0.1) at the cor¬ 
responding 6. And that of the red(dai'k yellow) empty circle denotes 
the average value of |m/| for Simp = 0.025(0.1). 

required to be unity for renormalization condition. In 
order to convergence of our calculations we have ver¬ 
ified that further sampling almost does not change pA 
and p|,„|. 

3. weak scattering potential cases 

For weak SP, we choose typical doping concentra¬ 
tions to investigate the distribution of the order param¬ 
eters. In the figures of distribution function, the hor¬ 
izontal axis denotes the value of the normalized order 
parameters Q-JQo, and the vertical axis denotes the cor¬ 
responding probability pQ, where Q stands for A or \m\. 

FigHJis for Uimp = 0.5 with 6 imp chosen as 0.025 
and 0.1. For 6 imp - 0.025, the peak of pA is located 
at Ai/Ao = 1.0, despite the value of 6 . It means that 
the most probable value of SC order is Aq. Black 
solid line in FigUfa) shows that for underdoped case 
6 - 0.04, the value of Ai/Ao is in the range [0.85,1.07]; 
for larger 6 = 0.08, a narrower range can be seen in 
FigHfb). Panel(c) plot the overdoped cases 6 - 0.16, in 
which the magnetic order vanishes and the peak of pA is 
sharper. Increase 6 imp to 0.1, the corresponding pA be¬ 
comes wider with lower peaks which are plotted by the 
blue dotted lines in FiglT] 

Magnetic order is more sensitive to multi-impurity, 
the peak of distribution function p|„,| is lower than that 
of corresponding pA. As Simp - 0.025, for S - 0.04 the 
range of |OTi|/|mo| is [0.95,1.15], for large doping level 
S - 0.08 it is widely extended to [0.8,1.2] as shown by 
the red dashed lines in Figina)(b). For Simp = 0.1, the 
corresponding pfmi is becomes wider than that of Simp - 
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Figure 2: (color online) (a)Plots of pA as a function of Aj/Ao for dif¬ 
ferent Simp, at a fixed set of (5 = 0.04, Uimp = 2.0. (b) Similar to 
(a) but for p|,„|. From right to left the black solid, red short dashed, 
blue short dotted, pink dash-dotted, and green dashed lines represent 
Si„,p = 0.0125,0.025,0.05,0.1,0.2, respectively. 

0.025 which are shown by the dark yellow dash-dotted 
lines in Fig[T] 

X coordinate of the square(circle) dots in Fig[T]is the 
average value of A,(|ot,|) in the corresponding case. As 
can be seen Q\ are close to Qq in all cases. Thus the 
system is very robust against weak SP regardless of 
Simp- Even when Simp - 0.2 (not shown here), Q\ does 
not have qualitative change with respect to Qo- From 
FigHtc), we note that for large S, pa has a small hump 
at the right side of Ai/Ao = 1.0. It is understandable 
since around the impurities Ai are enhanced for high S. 

4. moderate scattering potential cases 

Then we investigate moderate SP cases, and Fig]^ 
shows that for Uimp - 2.0, the corresponding pA and 
Pjmj extend much wider than that of Uimp = 0.5. As 
the impurity concentration Simp increases, the peak po¬ 
sition moves to smaller value and the peak height falls 
down as seen in FiglS^b) and all curves in Figl^Ja) ex¬ 
cept for the Simp = 0.2 case. For much big Simp, Qi on 
every site may be less than Qo, and the distribution will 
have a peak near zero with a considerable height. The 
green dashed curve in Figl^a) presents this regime with 
S/mp ~ 0.2. 

Distribution functions for S - 0.08 have similar prop¬ 
erties as that of 5 = 0.04. Figl^a) shows that as Simp 
increases gradually from zero to 0.2, the peak height 
of corresponding pA decreases. Since no one site has 
A; = 0, the system is in a stable SC phase for Simp - 
0.2, Uimp = 2.0. As for p|m|, the height of the peak 
goes up when Simp >0.1 which are shown in FigHb). 
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Figure 3: (color online) (a)Plots of pA as a function of Aj/Ao for dif¬ 
ferent 6imp, at a fixed set of 5 = 0.08, Uimp = 2.0. (b) Similar to 
panel a but for p\,„\. Inset of (b) plot p|„j| for Ni^p = 1. (c)Plots 
of Pa as a function of Aj/Ao for different Simp, at a fixed set of 
S = 0.16, Uimp = 2.0. From right to left the black solid, red short- 
dashed, blue short-dotted, pink dash-dotted, and green dashed lines 
represent Simp = 0.0125,0.025,0.05,0.1,0.2, respectively. 


(a) 5=0.04,6 =0.0125,U. =2.0 (b) 6=0.08,6. =0.0125,U. =2.0 
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Figure 4: (color online) Images of lwi|/|mol for Simp = 0.0125, Uimp = 
2.0. Panel (a) is for d = 0.04 and panel (b) is for S = 0.08. 

When pq enters this regime, the corresponding order Q 
tends to disappear. It is worth noting that different from 
d = 0.04, see the black solid line in Figl3lb), p|m| for 
d = 0.08 has two humps instead of one sharp peak. The 
reason is that when d near optimal doping, small amount 
of impurities can induce the lattice separate into two 
sublattice since each unit cell contains two inequiyalent 
Fe ions. This is consistent with previous study lll5n . For 
d = 0.08, p|m| of single impurity is plotted in the inset of 
panel(b) with two pronounced peaks shown up. Multi¬ 
impurity will smear out the two humps. In order to see 
it clearly, we show the images of normalized magnetic 
order |OTi|/|mo| for d = 0.04 and 0.08 in the following. 

We can see form Fig|4]that \mi\ at or around the impu¬ 
rities is enhanced for both doping concentrations. Panel 
(a) shows that for 6 - 0.04, the main part of the lattice 
has \m\\ less but close to |mo| leading to the sharp peak in 
FiglUb). While FiglUb) shows that for 6 = 0.08, except 
for the enhanced region of |ot,|, two sublattice crossed 
each other explicitly, which contribute to the two humps 


(b)6=0.08,6.„ =0.1,U.„ =2.0 
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Figure 5: (color online) Images of mj for Simp = = 2.0. 

Panel (a) is for S = 0.04 and panel (b) for S = 0.08. 



Figure 6: (color online) Images of Aj/Ao for Simp = 0.0125, Uimp = 
2.0. Panel (a) is for S = 0.04 and panel (b) for S = 0.16. Where the 
black solid squai'es stand for the sites of the multiple impurities. 

in FigOb) for 6 i,„p = 0.0125. 

Instead of |mi|/|mo|, Fig|5] depicts mi for 6 imp - 0.1, 
Uimp = 2.0 cases. The stripe-like pattern is still stable 
for d = 0.04 as shown in FiglSja). When the doping 
level increases to d = 0.08, nij = +0.04 for the pure 
sample; however, the multi-impurity makes the stripe¬ 
like pattern almost disappear, since we have many sites 
mi = 0 in FiglSjb). Whereas magnetic order is more 
sensitive to multi-impurity in comparison to no one site 
has zero SC order. 

Results of overdoped 6 = 0.16, moderate Uimp - 2.0 
are shown in Fig|3lc). The peak height of pA is much 
higher than that of d = 0.08 and S = 0.04 especially 
for lower impurity concentration. The main feature is 
very similar to that of d = 0.08 except that pA can have 
satellite peaks located at ^ > 1.0 fori5,„,p < 0.05. Thus 
in Fig|6]we present the image of Aj/Ao for d = 0.16 and 
6 - 0.04 for comparison. 

Figl6]shows that the multi-impurity effect on Aj is sig¬ 
nificantly different for low doping level and high dop¬ 
ing level. Panel (a) shows for underdoped 6 = 0.04, 
around the impurities A, is suppressed. Away from the 
impurities A, is close to Ao and corresponding to the 
sharp peak of Figl2ja). On the contrary, for overdop¬ 
ing level d = 0.16, on the impurity sites as well as 
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on most of their nn sites Ai are enhanced, which can 
be clearly seen in Fig|6tb). The sites of enhanced A, 
lead to additional satellite peaks located at > 1.0 
in FiglSjc). As dimp increases, interference between the 
multi-impurities will smear out these small peaks. 

Although order parameters may be enhanced at some 
sites, the average value Qi is less than the corresponding 
Qo for all cases. FiglTja) shows that A,- decreases as 
Simp increases for all doping levels. The decreasing rate 
depends on 6 as well as on Simp- All A,- have finite values 
at 6 imp = 0.2. Magnetic order |m,| is also reduced as 
Simp increases as shown in FiglTJb). For S - 0.08, |m,| 
has a small finite value at Simp =0.1, while the SDW 
pattern is already largely broken for Simp = 0.1 as seen 
in FiglJtb). 

FiglTJc) plots Aj as a function of T for different Simp 
at fixed Uimp - 2.0, and 5 = 0.04. The black dotted line 
denotes the impurity-free case for comparison. From 
top to down the curves for increased Simp have similar 
shapes. At low T, the curves are flat, and the higher 
Simp the shorter the flat region. As T increases. A, falls 
more quickly and vanishes at the critical temperature Tc- 
At lower Simp, Tc is less deviated from the clean limit; 
while for higher Simp, Tc is remarkably suppressed. 

5. strong scattering potential cases 

Strong SP multi-impurity suppress the order parame¬ 
ters more easily, the relevant properties are depicted in 
FiglDfor large Uimp - 5.0. We only present the results 
of ^ = 0.08, since other doping cases have similar re¬ 
sults. Panel(a) shows that when Simp ^ 0.05(blue short- 
dotted line), the order parameters are concentrated near 
zero with a very sharp peak for Simp - 0.1 (pink dash- 
dotted line), and the system totally loses superconduc¬ 
tivity. Magnetic order is even more sensitive to impu¬ 
rities, panel(b) shows that when Simp > 0.0125, sharp 


Figure 8: (color online) (a)Plots of as a function of Aj/Ao for 
different Simp, at a fixed set of d = 0.08, Uimp = 5.0. (b) Simi- 
lai' to panel(a) but for p|,„|. From right to left the black solid, red 
short-dashed, blue short-dotted, pink dash-dotted represent Simp = 
0.0125,0.025,0.05,0.1, respectively. Inset of panel(a) plots Aj as 
function of Simp for different S at Uimp = 5.0. Inset of panel(b) plots 
|m/| as function of Simp for different S at Uimp = 5.0. 

peaks of p\m\ are all squeezed to zero with vanished mag¬ 
netic order. 

For all doping levels, average order parameters de¬ 
crease more quickly as Simp increases which can be seen 
in the inset of Figl^a). It shows that at Simp = 0.1, we 
have A, = 0.0 for all cases and inset of Panel(b) shows 
that |m,| decreases sharply as well. Although at i5 = 0.04 
and Simp - 0.1, |ot,| has a very small finite value, the 
SDW pattern is broken. 

6. summary 

Motivated by experiments of TM doped iron-based 
superconductors and the different superfluid densities 
between bulk and film samples, we investigate multi¬ 
impurity effect systematically based on a phenomeno¬ 
logical model. s±-wave superconductivity has its dis¬ 
tinct properties when impurities are introduced. We use 
a number of configurations with random impurities to 
study the probability distributions of the order parame¬ 
ters. 

For weak SP, averaged order parameters do not have 
qualitative change even at very high impurity concen¬ 
tration Simp- They resemble the properties of the clean 
sample, just like by setting the Coulomb interaction U 
slightly larger and the system is still in a stable SC 
phase. 

For moderate SP, small Simp but larger S, the distri¬ 
bution functions of pA have a sharp peak together with 
a small hump, since the vicinity of impurities consists 


5 
























a considerably main part with enhanced A,. As 6 imp 
increases, the hump disappears due to the interference 
between many impurities. 

The multi-impurity effect on order parameters is re¬ 
markably different for low and high doping levels; how¬ 
ever, the averaged order parameters are suppressed with 
the increasing 6 imp in all cases. The probability distri¬ 
bution functions have similar behavior as dmp increases 
for moderate and strong SPs. The distribution peaks 
move towards zero for the larger 6 imp, meanwhile the 
peak near unity is suppressed and then replaced by the 
increased peak corresponding to a vanished order. 

s±-wave superconductor is robust again moderate SP 
multi-impurity effect, since A, is finite for all doping 
levels at large Simp - 0.2 and no site has A, = 0. While 
magnetic order is more sensitive to multi-impurities. At 
S - 0.1, although average |m;| is finite, the stripe-like 
SDW pattern has already been destroyed. The Tc val¬ 
ues of the samples are robust against small Simp, and are 
increasingly suppressed for the higher impurity concen¬ 
tration. 

Finally, we find that the strong SP is able to easily 
break down the order parameters. All the above features 
from our calculations provide some deep understand¬ 
ing of the difference between clean and dirty samples of 
iron-based superconductors. 
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